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ABSTRACT 

We report the discovery of an unambiguous, substantial low-energy cutoff in the broad band X-ray spectrum of 
the radio quasar RX J 1028. 6-0844 at a redshift of 4.276 obtained with the ASCA satellite, which we preferably 
explained as indication of excess X-ray absorption. The equivalent hydrogen column density of the absorbing mat- 
ter, depending on the redshift and metallicity, ranges from 2.5 10 21 cm -2 for local absorption up to 2.1 10 23 cm -2 
(solar metallicity) or 1.6 10 24 cm~ 2 (10% solar metallicity) for absorption at the quasar redshift. Such a value is 
among the largest found for high-redshift radio quasars. The absorption, if interpreted as being produced close to 
the quasar, may indicate the presence of a remarkably large amount of obscuring matter in the quasar environment 
in the early universe. Implications of the result for the possible origins of the absorbing matter are discussed, 
concerning especially galactic intervening matter, cool intracluster gas, and ambient medium around the quasar 
jet. The quasar itself has an enormous apparent luminosity of at least about 2.6 10 47 ergs" 1 {Ho = 50kms _1 Mpc" 1 , 
qo = 0.5) and a power law photon index of 1 .67(^ 04) in the 2-50 keV band in the source rest frame. 

Subject headings: X-rays: ISM -X-rays: galaxies - quasars: individual: RXJ1028. 6-0844 



1. INTRODUCTION 

Photoelectric absorption in the soft X-ray spectra of quasars 
provides unique probes of circum-source matter and quasar 
evolution. There is growing evidence that radio-loud quasars 
at high redshifts (z > 2) are commonly obscured by opaque 
matter — the X-ray spectra flatten toward low energies and this 
is attributed to absorption of the soft X-ray photons in excess of 
that due to Galactic interstellar medium, known as the excess 
absorption (Wilkes et al. 1992; Elvis et al. 1994; Serlemitsos 
et al. 1994; Siebert et al. 1996; Cappi et al. 1997; Reeves et 
al. 1997; Fiore et al. 1998). It is not known what cause such 
absorption, nor where they occur. Both, cosmologically in- 
tervening and quasar-associated materials are proposed as ab- 
sorber candidates; and the latter is favored by statistical argu- 
! ments (Bechtold et al. 1994; O'Flaherty and Jakobsen 1997; 
Fiore et al. 1998; Yuan and Brinkmann 1999) for high-z radio- 
quiet quasars showing no excess absorption, and by tentative 
indications for variability of the absorption (e.g. Schartel et al. 
1997, Cappi et al. 1997). This picture has been established 
based mainly on observations of quasars in the redshift range 
of 2-A. The highest redshifts up to which the excess absorp- 
tion was reported in these observations are of PKS 2126-158 at 
z = 3.3 (Elvis et al. 1994) and of S5 0014+81 at z = 3.4 (Cappi 
et al. 1997), and probably of Q1745+624 at z = 3.9 (Kubo et al. 
1997). 

At even higher redshifts (z > 4), only two quasars, both 
radio-loud, were observed with fair broad band X-ray spec- 
troscopy, 1508+5714 at z=4.3 (Moran and Helfand 1997) and 
GB 1428+4217 at z=4.7 (Fabian et al. 1998). While the 
ASCA spectra of both objects do not require excess absorp- 
tion, GB 1428+4217 was recently reported to have excess ab- 
sorption with a ROSAT PSPC observation (a hydrogen column 



density of 1.52(±0.28) 10 22 cm" 2 , assuming the absorber asso- 
ciated with the quasar; Boiler et al. 2000). Another interest- 
ing and probably related aspect is that both objects show blazar 
properties, i.e. radio and X-ray emission dominated by relativis- 
tically beamed components from the jets. Though two objects 
are too few for any statistical inference, it is intriguing to won- 
der whether this behavior is typical at very high redshifts or 
merely selection effects. 

In this paper, we report an ASCA observation of the z > 4 
quasar RXJ1028. 6-0844, as part of our X-ray high-z quasar 
program (Matsuoka et al. 1999). The quasar shows an unam- 
biguous, strong low-energy cutoff in its X-ray spectrum. The 
X-ray source RX J 1028. 6-0844 was discovered in the ROSAT 
All-sky Survey (RASS) and identified with a quasar at z = 4.276 
by Zickgraf et al. (1997). Its X-ray colors in the ROSAT en- 
ergy band indicated an extremely hard spectrum, thus, strong 
X-ray absorption if the 'intrinsic' spectrum is of typical radio- 
loud quasars (Zickgraf et al. 1997). The object is also associ- 
ated with a radio source, PKS B 1026-084 (Otrupcek & Wright 
1991), which has a flux density of 220 mJy at 5 GHz and a flat 
radio spectrum (a = 0.3, S oc v~ a ); hence it is a flat-spectrum 
radio-loud quasar (FSRQ), and has an intense 5 GHz radio lu- 
minosity 5 vL v = 5 10 44 ergs _1 . We present the observation and 
data reduction in Sect. 2, and spectral and temporal analyses in 
Sect. 3 and 4, respectively. In Sect. 5 the implications for the 
origins of the excess absorption are discussed. A summary of 
the main results is given in Sect. 6. Errors quoted are of the 
68% confidence level (c.l.) throughout the paper unless men- 
tioned otherwise. 



2. OBSERVATION AND DATA REDUCTION 
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RXJ1028.6-0844 was observed with the ASCA satellite 
(Tanaka et al. 1994) on Nov. 25, 1999 and the duration lasted 67 
hours. The Solid-state Imaging Spectrometers (SIS) were oper- 
ated in the 1-CCD faint mode and the Gas Imaging Spectrom- 
eters (GIS) in pulse height mode. For SIS the 'bright 2' mode 
data (grade 0, 2, 3, 4) were used, which were converted from 
the original faint mode data. Corrections for the dark frame er- 
ror and the echo effect were applied. The data reduction was 
performed in the standard way by using the FTOOLS (v.4.2) 
utilities. Hot and flickering pixels were removed from the SIS 
data. To improve the reliability of the results, we chose conser- 
vative data screening criteria, which are listed in Table 1 . 

The source counts were extracted from circular regions of 
3.5 and 5 arcmin radii for SIS and GIS, respectively. The back- 
grounds were determined in two ways and compared with each 
other: from the blank-sky observations using the same region of 
the detector, and from source-free regions in the field of view 
of the same observation (local). For the latter, for each of the 
GIS detectors two local off-source regions were used to extract 
background events, which were chosen so as to have the same 
radii and off-axis angles as the source region; for SIS, the back- 
ground events were extracted from the whole CCD chip with a 
circular region of 5 arcmin radius around the source excluded. 
In addition to closely matching the screening criteria, the lo- 
cal backgrounds have the advantages over the blank-sky back- 
grounds regarding the uncorrected decreasing of the low-energy 
efficiency for SIS and the variability of the internal background 
for GIS. This is especially true considering the long time span 
between the blank-sky and the current observations. We thus 
mainly used the local backgrounds to derive results and used 
the blank-sky backgrounds for comparisons. The resulting ef- 
fective exposures and extracted source counts are given in Ta- 
ble 1. 

Since the two GIS have identical response matrices and the 
two source spectra agree well with each other, we combined 
them into one single GIS spectrum 6 . All the spectra were re- 
binned to have at least 30 counts in each energy bin. We used 
the most updated releases of the response matrices available so 
far (V4.0 GIS RMF gis2(3)v4_0.rmf, VI. 1 SIS RMF and the 
calibration file sisph2pi_l 10397. fits) in spectral fit. 

3. SPECTRAL ANALYSIS 

It is known that the responses of the two SIS suffer from cali- 
bration uncertainties, which are the most problematic at around 
0.6 keV (ASCA Guest Observer Facility web page). To min- 
imize the possible systematic error introduced in the results, 
we ignored the SIS data below 0.7 keV in spectral fit, but com- 
pared them with the fitted models. The energy ranges used for 
spectral fits are 0.7-9.5 keV for both SIS and GIS. We first fit 
the spectra with a simple power law model modified by photo- 
electric absorption assuming an absorber in the observer frame 
(local absorption, Sect. 3.1) and at a redshift (redshifted absorp- 
tion, Sect. 3.2), respectively, in addition to the Galactic absorp- 
tion. Then we test other spectral models for the continuum in 
Sect. 3.3. We consider the impact of the excluded low-energy 
data on the obtained results in Sect. 3.4. 



gle power law model absorbed by 'cold' gas (model wabs 1 in 
XSPEC v.10) with the column density fixed at A^ al gave a flat 
photon index T = 1.34 ± 0.03 (joint fit of the SIS0, SIS1, and 
GIS spectra, see Table 2 for the fitted parameters). Though 
the model is statistically acceptable, the data below 1 keV fall 
systematically below the model. The deviations are more pro- 
nounced as shown in Fig. 1, in which the above model was fit- 
ted to the spectra in the restricted 2.5-9.5 keV band only (13- 
50 keV in the quasar frame) and extrapolated down to 0.5 keV. 
Such a low-energy cutoff cannot be explained by the known 
calibration uncertainty for SIS. This is because, firstly, the GIS 
spectra, though being less sensitive in the low energy range, 
show the same effect with similar amounts of deviations as the 
SIS spectra; secondly, the loss of the SIS efficiencies is quanti- 
fied as up to about 20-40% at 0.6 keV (ASCA Guest Observer 
Facility web page), while in this case the deviations of the data 
from the model are about a factor of two on average, for both 
SIS and GIS. 

Setting the absorption column density A^h as a free parame- 
ter improved the fits significantly, and yielded large Nh for both 
SIS and GIS. The x 2 value was reduced by A x 2 =6, 18, and 34 
for fitting the GIS, SIS0+SIS1, and the joint GIS+SIS spectra, 
respectively. Thus, the addition of the excess absorption term is 
significant at > 99.99% c.l. For the joint GIS+SIS fit this gave 
Nft = 3.0£%p 10 21 cm" 2 , and a photon index T = 1.67C$g), 
which is typical of FSRQ in the ^1-10(20) keV band at low and 
medium redshifts (e.g. Siebert et al. 1996; Lawson and Turner 
1997; Reeves et al. 1997; Cappi et al. 1997; Brinkmann et al. 
1997). The residuals of the fit are shown in Fig. 2. The SIS 0.5- 
0.7 keV data, which were excluded from the spectral fitting, are 
also plotted; it is shown that the extrapolation of the model does 
not under-predict the count rates in this band (though the real 
count rates could be somewhat higher because of the loss of 
the SIS efficiencies). Fig. 3 shows the confidence contours for 
two interesting parameters for the fitted A^h and T (marked by 
z = 0). Absorption of the quasar X-rays in excess of that due to 
the Galactic column density is evident. Using the background 
spectra derived from the blank-sky observations gave consistent 
results. 

Due to the increasing degradation of the response, the cal- 
ibration for both SIS is uncertain below ~ 1 keV; this effect 
reportedly results in deficits of the observed counts at low en- 
ergies and, effectively, apparent 'excess absorption' . From the 
most recent quantification of this problem (Yaqoob et al. 2000), 
the apparent excess Nr at the time of the observation was esti- 
mated as (using the 0.5-10keV data) ANn ~ 0.7 10 21 cm" 2 for 
SIS0 and about 1.0 10 21 cm" 2 for SIS1. Fitting the whole 0.5- 



lOkeV SIS spectra yielded JV H = 2.9(!^) 10 



,21 



3.1. Power law with local absorption model 

The Galactic column density in the source direction is A^ al 
= 4.59 10 20 cm" 2 (Dickey & Lockman 1990). A fit with a sin- 

6 The counting errors in each bin were assigned as pure Poissonian errors of the summed counts. 

7 For which the photoelectric absorption cross sections from Morrison and McCammon (1983) are used. 



cm" 2 (r= 1.62) 
for SIS0 and 4.5(!^) 10 21 cm" 2 (T= 1.74)for SIS 1; the result- 
ing excess absorption column densities are significantly larger 
than the estimated systematic uncertainty. Thus, the observed 
low-energy cutoff is a true feature of the source spectrum, not 
an artifact due to the calibration problem. Very roughly, we 
estimated the 'true' absorption column density by subtracting 
the systematic excess ANh from the fitted value, as Nu~ANh- 
We found that, for both SIS, the remaining column densities 
are very close to the fitted A^h values by using the 0.7-10keV 
spectra, for which the calibration uncertainty is thought to 
have less effect. Independently, the GIS spectrum gave a Nh 
3-OCjg) 10 21 cm" 2 in good agreement with that fitted using the 
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0.7-10 keV SIS spectra (see Table 2), suggesting that the SIS 
calibration systematic uncertainties are small and negligible in 
the 0.7-10keV band in comparison to the statistical uncertain- 
ties. We therefore regard the spectral parameters derived by 
using the 0.7-10keV SIS spectra, as performed throughout this 
work, as being little affected by the calibration uncertainties, 
and ignore this effect in the following analyses. 

For the power law with Galactic absorption model, the ab- 
sorption corrected 2-10keV flux is 1.2 10~ 12 ergs" 1 cm" 2 , cor- 
responding to a luminosity of 6.2 10 46 (2.4 10 47 )ergs _1 in the 
2-10 (2-50) keV band in the quasar rest frame. For the free 
absorption model, the luminosity corrected for the total (Galac- 
tic + excess) absorption is 9.8 10 46 (2.6 10 47 )ergs _1 in the 2- 
10 (2-50) keV band. These values were obtained using the GIS 
data, which are slightly higher than the SIS measurements by 
less than 10%. No iron K emission line was detected; the upper 
limit of the equivalent width for a narrow (a = 0.01 keV) line at 
6.4 keV is 195 eV at the 90% c.l. in the source rest frame. We 
found no iron K-shell absorption edge feature for local absorp- 
tion, setting a 90% upper limit on the optical depth as r < 0.79 
at7.1keV. 

3.2. Power law with redshifted absorption 
3.2.1. Absorption at the quasar redshift 

Assuming the absorber is associated with the quasar, we fit- 
ted the spectra with models of a power law modified by absorp- 
tion occurring at z = 4.276, and further by the (fixed) Galac- 
tic absorption. For 'cold' absorption models this results in the 
same \ 2 with, but a larger excess absorption column density 
A^ xc (2.1 10 23 cm~ 2 , see Table 2) by two orders of magnitude 
than, those for the above local absorption model. This is be- 
cause the absorbed photons have increased energies (by a fac- 
tor of 1+z) and thus decreased absorption cross sections in the 
absorber frame. We over-plotted in Fig. 3 the confidence con- 
tours for the fitted Ng c and T (marked by z = 4.28). It should 
be noted that we have assumed the solar elemental abundances 
for the absorbing gas. The energy range of E > 0.7 keV in 
which the X-ray photons are absorbed in the ASCA bandpass 
corresponds to E > 3.5 keV in the quasar rest frame; at such 
energies the cross section of photoelectric absorption is dom- 
inated by heavy elements, mainly O, Ne, Si, S, Ca, and Fe. 
Since at high redshifts the metal abundances could be well sub- 
solar (about 10% or less, e.g. Lu et al. 1996; Pettini et al. 1997; 
Prochaska & Wolfe 2000), we assumed the metal abundances 
as A z =0.1 A ZQ and repeated the fit (using zvfeabs model in 
XSPEC). The resulting absorption column density becomes as 
high as Afjf° ~ 1 .6 10 24 cm" 2 (Table 2). Since the true elemental 
abundances are unknown, the determination of N^ c is subject 
to a large uncertainty, however. 

We searched for any redshifted iron K absorption edge by fit- 
ting the edge separately from the redshifted absorption model. 
The best-fit edge energy is 7.2(^ f ) keV; the optical depth at the 
threshold energy is t = 0.15(^ J 2 ,), which is consistent the pre- 
diction from the fitted A^ xc (r = 0.26 for A Z =A Z0 ). However, 
given the reduction in \ 2 of A\ 2 = 1.3 only, the detection of 
the iron K edge cannot be justified. The non-detection of a sig- 
nificant iron edge gave no determination of the redshift of the 
absorber, unfortunately. 

It is possible that the absorber is photoionized by the intense 
radiation from the quasar. A photoionized absorption model 
(absori in XSPEC) was fitted to the spectrum, in which the 
power law X-ray emission of the quasar was taken as the ion- 



izing continuum. The result is inconclusive, however. The fit 
converged at neutral absorption with the best-fit ionization pa- 
rameter (as defined in Done et al. 1992) £ ~ 0, while ionization 
of the gas to some degree is allowed, £ < 2 (20) at the 90% 
(99%) c.l. For sub-solar metal abundances the allowed £ range 
becomes lower, £ < 0.1 at the 90% c.l. for A z =0.1A ZQ . It 
should be noted that, however, if the relative abundance of iron 
with respect to other heavy elements is much lower than the so- 
lar value, the ionization state is almost unconstrained, though 
this seems to be an unlikely case. For instance, for A z =A ZQ 
except A Fe = 0.1 A Fe 0, the 90% confidence range for £ could be 
as large as 0-1000. This is because the constraint on the ioniza- 
tion state comes largely from the lack of the predicted K-edge 
features of highly ionized irons in 1.3-1. 7 keV in the observer 
frame. 

3.2.2. Redshift dependence 

The unknown location of the absorber introduces a large un- 
certainty in the determination of A^ xc . With the limited low- 
energy bandpass and energy resolution, ASCA cannot distin- 
guish between local and highly redshifted absorption models. 
This is because, firstly, these models could have similar spec- 
tral shapes of the absorbed continuum in the ASCA band, and 
secondly, the photon statistics are generally poor at low ener- 
gies anyway, especially in the case of strong absorption. Leav- 
ing the redshift free yielded no statistically significant range in 
the parameter space, i.e. similar goodness-of-fits for redshifts 
in the full range of 0^.28, though there exist a few shallow 
local minima in the \ 2 distribution (however, see Sect 3.4 for 
the case when the 0.5-0.7 keV spectral data are included). The 
dependence of the excess column density on the assumed red- 
shift is shown in Fig. 4, in which the confidence contours are 
plotted on the -z plane. The required A^f c decreases down 
to 2.5 10 21 cm -2 for z=0, as shown in Sect. 3.1. It is also shown 
that, at a given redshift, A^ xc is relatively well constrained. It 
should be noted that these results were obtained by assuming 
the solar elemental abundances; if the metallicity is signifi- 
cantly sub-solar, the inferred A^f c value would increase sub- 
stantially (see above). 

3.3. Other continuum models 

Thermal bremsstrahlung model — We also fitted the spec- 
trum with redshifted thermal bremsstrahlung plus local ab- 
sorption models. For absorption fixed at the Galactic A^ al , a 
source frame temperature as high as 1 32(^7) keV is required. 
Free fitting the absorption column density gave T = 55Cg 2 ) keV 
and, similarly, a large Nn =2.1Cq 4 ) 10 21 cm" 2 . The resulting 
goodness-of-fit is the same as that for a power law fit. For such 
a high temperature, the thermal bremsstrahlung spectrum is in- 
distinguishable from a power law one in the ASCA band. 

Broken power law — A convex spectral shape, like what ob- 
served in some BL Lac objects, can result in apparent excess 
absorption when fitted with a power law model. Such spectra 
are usually parameterized with a broken power law model. Fit- 
ting a broken power law with fixed Galactic absorption yielded 
statistically acceptable fit, x 2 = 190 for 209 degrees of freedom 
(d.o.f.); in spite of one more free parameter, the \ 2 is some- 
what worse than those for the power law with excess absorp- 
tion models. The resulting break energy of 9 ± 1 keV (in the 
quasar frame) is higher than those found for the X-ray spectra 
of BL Lac objects (a few keV or less, e.g. Barr et al. 1988; 
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Madejski et al. 1991). Leaving the absorption Nu free im- 
proved the fits slightly (A% 2 = 3), and yielded similar values 
for the high- and low-energy photon indices (Table 2), and a 
N H ~ 2.8(±0.7) 10 21 cm" 2 significantly larger than the Galactic 
value. Such a model is effectively consistent, within the uncer- 
tainties of the fitted parameters, with the absorbed power law 
models — the best-fit models as shown above. 

3.4. Effect of the low-energy end data: excess absorption? 

As shown above, the 0.7-9.5 keV spectrum can be similarly 
well modeled by either the power law with (local or redshifted) 
excess absorption models or the broken power law continuum 
with Galactic absorption. Since these models show increas- 
ing spectral divergence toward low energies, the omitted 0.5- 
0.7 keV SIS data, even in the presence of the calibration uncer- 
tainties, might be useful for discriminating between these mod- 
els. Here we examine the effect of the excluded 0.5-0.7 keV 
SIS data (two bins after spectral binning) on the obtained re- 
sults in two ways. 

Firstly, we extrapolated the above best fits for various mod- 
els obtained with the 0.7-9.5 keV spectrum down to 0.5 ke V for 
SIS and compared the model values with the data. For all the 
models, deficits of the observed count rate in these energy bins 
were always found for both SIS0 and SIS 1, which might be a 
result of the above systematic effect, or mismatched models, 
or both. We examined the consistency of the observed 0.5- 
0.7 keV count rates with the current models by evaluating the 
F-statistic F, which basically measures the change in x 2 with 
the change of the degrees of freedom in fitting. (For a F as 
large as exceeding the F a which corresponds to a certain prob- 
ability level P(F CI ), the chance probability under the model hy- 
pothesis is as small as P < P(F cr )). For models of the power 
law with redshifted absorption, with local absorption, and the 
broken power law with Galactic absorption, we found the \ 2 
increments of A% 2 = 4, 12, and 30 for two additional d.o.f., 
and the calculated F=2.l, 6.4, 14.4, and the probability levels 
P > 5%, < 1%, < 1%, respectively. Thus, regardless of the SIS 
calibration uncertainty, the second and the third models do not 
agree with the 0.5-0.7 keV data at high confidence levels. This 
can be understood as the redshifted absorption model, with its 
large absorption Nh, predicts the strongest low-energy spectral 
cutoff, i.e. the least count rates in 0.5-0.7 keV. 

It has been claimed that at 0.6 keV the decrease of the effi- 
ciencies can be as much as 40% for SIS 1 in the current observa- 
tions (compared with GIS, ASCA Guest Observer Facility web 
page). We took this systematic effect into account by lowering 
the model values in the 0.5-0.7 keV bin by 40% for (conser- 
vatively) both SIS0 and SIS 1, and repeated the above testing. 
We have A% 2 =1, 5, and 9.5 for the above three corresponding 
models, respectively. This means while the first and the second 
models give the best and marginal consistencies with the 0.5- 
0.7 keV data (F=0.6 and 2.7), respectively, the broken power 
law model still disagrees with the included data at high confi- 
dence level (F=5, P < 1 %). 

Secondly, we tried to compensate the SIS calibration un- 
certainties by decreasing the detecting efficiencies in the 0.5- 
0.7 keV band by 20% for SIS0 and by 40% for SIS1, and fit- 
ted the SIS+GIS spectra in the whole 0.5-9.5 keV band using 
the 'corrected' SIS calibration files. This resulted in the best- 
fit model being the power law with excess absorption in the 
quasar frame, with the goodness-of-fit (x 2 =l90 for 212 d.o.f.) 
improved slightly over the power law with local absorption 



model (x 2 =193 for 212 d.o.f.), and considerably over the bro- 
ken power law with Galactic absorption model (x 2 =197 for 21 1 
d.o.f.). The fitted spectral parameters are similar to those ob- 
tained by fitting the 0.7-9.5 keV spectra. It is interesting to 
note that, for the excess absorption models, the redshift of the 
absorber could be free fitted as z = 4.2^ 5 (90% c.l.); however, 
this result is only suggestive given the non-uniform efficiency 
loss over the 0.5-0.7 keV band and the less strict treatment here. 

We conclude that the broken power law with Galactic absorp- 
tion model gives a poor description to the observed low-energy 
spectrum down to 0.5 keV, in comparison with the (redshifted) 
excess absorption model. Furthermore, the SIS data in the 0.5- 
0.7 ke V range are better described by the absorption model with 
the absorber being at or close to the quasar redshift than the lo- 
cal absorption model. 

4. TEMPORAL PROPERTIES 

The 67-hour duration of the observation corresponds to ~ 
12.7 hours in the quasar frame, which is short for quasar tim- 
ing analysis. The X-ray variability was searched for from the 
background-subtracted light curve by using the packages pro- 
vided in the FTOOLS. No statistically significant variations 
were found in both the co-added SIS and GIS light curves, as 
well as in the combined SIS+GIS light curve. It is noted, how- 
ever, that a marginal, small amplitude variation could probably 
be present in the first half of the duration (x 2 test, % 2 = 15 for 9 
d.o.f., a probability level of 0.09). 

We compared our results with the RASS flux by extrapolat- 
ing the best-fit model down to 0. 1 keV. The predicted count rate 
in the ROSAT 0. 1-2.4 keV band is 0.02 ctss" 1 , consistent with 
the RASS measurement 0.035 ± 0.01 1 cts s -1 (Zickgraf et al. 
1997) within the errors. Thus, no significant flux variation was 
found by comparing the two observations of about 2 years apart 
in the source rest frame. It should be noted that, in the ROSAT 
soft 0. 1-0.41 keV band, only ^0.1 source count is predicted for 
the 424 second RASS exposure due to the strong absorption of 
soft X-rays, in contrast to the observed ~ 2 counts (a counting 
statistical probability of 3% only). This may imply the latter 
to be either of background fluctuations, or source photons scat- 
tered off the absorbing medium down to the soft X-ray energies. 

5. DISCUSSION 

We have shown that the X-ray spectrum of RX J1028. 6-0844 
exhibits a low-energy cutoff much strong than that expected 
from the Galactic absorption. Though both, the power law with 
excess absorption and the broken power law with Galactic ab- 
sorption models give acceptable fits to the spectra, the former 
is favored over the latter based on the fitting statistics. On the 
other hand, the power law model has been found to be a char- 
acteristic of the intrinsic X-ray continuum emission of FSRQ 
over a wide range of redshift. We thus suggest that the broken 
power law model seems to be an unlikely case, though it cannot 
be ruled out confidently, and the observed excess spectral cut- 
off is most likely caused by absorption of the soft X-ray photons 
from the quasar. This explanation is much natural in the light 
of what has been established in the redshift range of around 2 
to 3 — excess X-ray absorption is common in radio quasars, and 
extends this behavior up to redshift above 4 (two out of the three 
objects with well determined X-ray spectra). We consider this 
scenario only and discuss its implications in below. 

5.1. The absorber 
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The derived column density of the excess absorption (here- 
after denoted simply as Nu) is among the largest 8 ever found for 
high-z radio quasars, which imposes a strong observational con- 
straint on the nature of the absorbing material. In addition, any 
satisfactory absorption models should be able to accommodate 
the absorber to the following facts: (a) the relatively well con- 
strained Nu -z relation (Fig 4), though with some uncertainties 
introduced from the unknown metallicity; (b) being nearly neu- 
tral or moderately ionized, £ < 2 at the 90% c.l. (unless the iron 
abundance relative to other heavy elements is abnormally sub- 
solar; see above); (c) no strong absorption beyond the Lyman 
limit, as indicated by both the optical (UV in the source frame) 
spectrum and the measured B-magnitude (emitted at 834 A) in 
Zickgraf et al. (1997); and similarly, (d) no indication of heavy 
dust extinction; (e) if the X-ray emission is mainly from the jet 
(a likely case, see Sect. 5.2), the absorber must lie outside the 
jet from the nucleus or around the jet. Several possible origins 
of the X-ray absorption of high-z radio quasars have been pro- 
posed (see Elvis et al. 1994, 1998), such as damped Lya sys- 
tems, intracluster gas, and the boundary layer between the jet 
and surrounding medium. We discuss in below the implications 
of our results for the possible absorbers of both intervening and 
intrinsic origins. 

5.1.1. Cosmological intervening material? 

An inspection of the optical image of the quasar field ob- 
tained by Zickgraf et al. (1997) suggests readily an absorber 
candidate: an extended, faint (speculative) spiral galaxy, which 
lies ~ 7" away from the quasar and seems not highly inclined 
given the estimated optical axes of ~ 5" x 3". The presumed 
H I gas disk/halo might intersect the quasar line-of-sight at a 
transverse distance of from a few kpc (z <C 1) to several tens 
kpc (z ~ 1). If this object is indeed a spiral galaxy, absorp- 
tion of the quasar X-rays by its H I gas is unavoidable. The 
question is whether the amount of the absorbing gas can ac- 
count for the large column density derived here. Observational 
evidence for X-ray absorption due to galactic H I disk/shell 
has just emerged in a few systems (Elvis et al. 1997; Yuan 
and Brinkmann 1999), but the measured A^h are relatively low 
(several times 10 20 cm~ 2 ). We tried to assess quantitatively 
the amount of such possible absorption. We estimated the 
averaged A^h of the anticipated H I gas in this galaxy from 
its optical brightness, by making use of the galactic H I - 
mass and luminosity relation and the distance-independent scal- 
ing N H ~ 10 20 (M ra /LB)10 a4(27 -' tB) (Disney and Banks 1997), 
where Mm/L B is the H I -mass to light ratio in units of M Q /L Q , 
and /iB the Z?-band surface brightness averaged over the H I 
disk. Adopting the measured nig = 20.2 ± 0.5, we found a con- 
servative fi B ~ 25.2 (Galactic extinction corrected), which is 
averaged over the solid angle confined by an ellipse with the 
inferred axis ratio and intersecting the quasar sight-line. This 
led to an averaged Nu ~ 5.2(+] \) 10 20 (M H i/Lb) cm" 2 , where the 
quoted errors are from the uncertainty of m B . To match the ex- 
cess absorption N R of 2.5C$%) 10 21 cm" 2 (z = and A z = A Z0 ), 
a Mhi/Lb ratio as high as 4-5 is needed; this seems to be un- 
likely as spirals being found to have Mm/L B < 1 regardless of 
the morphological type (e.g. Kamphuis et al. 1996). Moreover, 
the true Nh at the site of the quasar sight-line might be even 
lower than the averaged value as the H I surface density drops 
rapidly outward from the center. We therefore suggest that this 
galaxy is probably responsible for part of the X-ray absorption, 

8 Similar amounts of excess absorption were also reported in S5 0014+81 (Cappi 



but it is unlikely — though we cannot rule it out — to be able to 
account for the total absorption opacity. After subtracting the 
plausible contribution of this galaxy (assuming Mm/L B < 1), 
an excess column density of Nu > 2. 0(^4) 10 21 cm" 2 (z = 
and Az = Az©) may remain unaccounted for. 

Such a large amount of absorption seems not to be easily ex- 
plained by the general damped Lya systems (DLA), neutral hy- 
drogen reservoirs with the largest column density (Nu i) known 
in the universe. Although DLA may have the distribution of 
A^h 1 reaching a few times 10 21 cm" 2 at its high end, the steep 
distribution function with a power law index of -3 renders such 
systems to be extremely rare (Wolfe et al. 1995, Zwaan et al. 
1999). The deficit of column density is even severe regarding 
the substantially sub-solar metal abundances (Az < O.IAz©) of 
DLA, even at z ~ (Miller et al. 1999), which require Nu to 
be well above 10 22 cm" 2 in order to produce the observed ab- 
sorption. Although the true total hydrogen column density A^h 
could be somewhat higher than A^h i considering a possible ion- 
ization fraction of H, the correction should not be large as gases 
in DLA are believed to be mainly neutral (Wolfe 1993). A de- 
tailed study by O' Flaherty and Jakobsen (1997) incorporating 
absorber statistics confirmed the extremely low detectability of 
high-A^H intervening systems; the chance probability to have 
a DLA with Nu > 10 21 cm" 2 on the line-of-sight to a high-z 
quasar is less than a few percent. This probability is even lower 
considering that any DLA systems, if do produce the absorp- 
tion, are most likely at low redshifts so as to avoid the otherwise 
much higher X-ray A^h required (Fig. 4). The non-detection of a 
DLA in the optical spectrum obtained by Zickgraf et al. (1997) 
places an upper limit on the redshift of any possible DLA as 
~ 2.5. Therefore, the likelihood of the absorption being caused 
by cosmologically intervening material is regarded to be low. 

5.1.2. Intrinsic absorber? 

Alternatively, the excess absorption might take place intrin- 
sically close to the quasar. This seems to be plausible in the 
light of the statistical arguments for high-z radio-quiet quasars 
showing no excess absorption. The inferred absorption column 
density in the quasar rest frame then reaches the highest ever 
found among this type, ~ 2 10 23 cm" 2 (> 1.7 10 23 cm" 2 cor- 
rected for the possible foreground galactic absorption assuming 
the expected Mui/L B < 1, see Sect. 5.1.1) or even higher up to 
~ 10 24 cm" 2 if Az < O.IAz©. Such thick gas can not be easily 
explained by the known mass components. This result strength- 
ens the tentative, apparent trend of the increasing Nu of excess 
absorption in radio quasars with the increase of redshift, which 
might have interesting consequences for understanding the evo- 
lution of quasars and, probably, the formation of the cosmic 
structure. 

One immediate problem, which seems to be common among 
objects of this type (Elvis et al. 1994), arises as the apparent in- 
consistency between the X-ray and optical data: the thick gas, 
if contains dust similar to that of our Galaxy in both composi- 
tion and content, would have blocked almost all the optical/UV 
light of the quasar believed to originate from the central re- 
gion, e.g. A(A1230A)=40 (Seaton 1979) for the observed R- 
magnitude. Though this problem may be relieved by a likely 
extremely low dust content or different dust compositions at 
high redshift (see Elvis et al. 1994 for a discussion), the lack 
of strong Lyman limit absorption remains a concern. A rough 
estimate of T 9n p^ ^ 1 from the optical spectrum in Zickgraf et 

et al. 1997) and PKS 0528+134 (Reeves et al. 1997). 
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al. (1997) suggests the absorber to be highly ionized, Nu i /Mi 
< 10" 6 (taking the conservative value of Mi =1-7 10 23 cm" 2 ). 
This ionization state requires the dimensionless ionization pa- 
rameter U > 5 (U is defined as the hydrogen-ionizing photon- 
to-electron number density ratio) for the typical AGN ionizing 
continuum, or, when converted to £, £ > 250; assuming the ion- 
izing continuum as the extrapolation of the observed flat power 
law with r ~ 1 .7 requires £ > 370. However, we found £ < 2 
(20) at the 90% (99%) c.l. from the X-ray data. One natural 
explanation would be that the UV continuum emitting region is 
only partially or not covered by the X-ray absorber. If this is the 
case, special geometry of the absorbing gases may be required; 
for instance, the gas might be ambiently local to the jet if the 
X-rays are of jet emission (constraint e). Again, it should be 
noted that, if the gas has iron-poor metallicity in comparison to 
the solar one, a higher ionization state cannot be ruled out. 

Interestingly, Elvis et al. (1994) related the X-ray absorption 
of quasars up to z ~ 3 to cooling flows of the possible host clus- 
ters of galaxies, though the existence of relaxed clusters at such 
high redshifts is unknown. If we take this hypothesis and ex- 
tend the redshift up to z - 4 for the case of RX J 1028. 6-0844, 
its large Nu value in comparison to those found from X-ray 
absorption in low-redshift clusters (~ 10 21 cm" 2 , White et al. 
1991) would imply an extremely strong cooling flow. The mass 
of the cooled gas reaches ~ 3 10 13 /? 2 00 M Q for a conservative 
Mi=10 23 cm -2 , where /?ioo is the cooling flow radius in units 
of 100 kpc. For R\qq ~ 1 as typically found in low-z clusters, 
this means a cooled mass of ~ 3 10 13 M Q , and a mass cooling 
rate > 3 10 4 M Q yr" 1 given the short life-time of the system of 
less than 10 9 years. Such a cooling flow almost approaches the 
'maximal cooling flow' in which the cooling time is compara- 
ble with the gravitational free-fall time (Fabian 1994). Though 
these numbers will be somewhat reduced for a likely clumpy 
distribution of the cooled gas, a cooling flow stronger than those 
in nearby clusters is suggested. 

It has also been noted that the massive (<~ 10 14 M Q ) neu- 
tral hydrogen cloud at z=3.4 reported by Uson et al. (1991) 
could be a possible origin of X-ray absorption. The gas cloud 
was detected via both absorption and emission lines at the ra- 
dio 21 cm wavelength, and was interpreted as indications of 
a proto-cluster of galaxies. The H I column density derived 
from the absorption line is 4.4 10 18 (r s /K) cm" 2 (T s the spin 
temperature of H I); the mean surface mass density was es- 
timated from the emission line as ~ 48 M pc~ 2 (Uson et al. 
1991), which corresponds to an averaged column density A^h i 
~ 5.7 10 21 cm" 2 . If such a system is responsible for the X- 
ray absorption in RX J 1028. 6-0844, the gas must be hot with 
T s > 10 4 K, and/or somewhat ionized, Nh /Mi i > 10; these 
lower limits will be further raised up for Az <Az©. Thus, the 
inferred T s and/or total gas mass are likely exceeding 10 4 K and 
3 10 14 M Q , respectively, the values proposed for a "Zel'dovich 
pancake" (Sunyaev and Zel'dovich 1975). 

If the quasar X-rays are dominated by jet emission, the ab- 
sorption may be caused by intimately jet-linked material. Such 
a scenario explains naturally the lack of UV absorption and 
reddening, which seems to be a commonplace for high-z ra- 
dio quasars showing excess X-ray absorption. It is of particu- 
lar interest to relate the X-ray absorption with the optical-radio 
'alignment effect' found in high-z radio galaxies, as noted in 
Elvis et al. (1998). This effect is the alignment of the ex- 
tended, polarized optical emission with the radio axes (e.g. 
McCarthy 1993), which is thought to be of (dust) scattering 



and/or star-formation origin. A class of dust scattering mod- 
els advocates the existence of gas-dust clouds embedded in a 
multi-phase boundary layer formed (from interaction) between 
the jet and ambient medium (e.g. De Young 1998). Although 
such a jet environment is largely uncertain, it is interesting to 
note the following order-of-magnitude estimation. We consider 
a specific set of model parameters which was adopted so as to 
roughly satisfy the physical condition and observations as pro- 
posed in De Young (1998): gas density n=10-100cm" 3 , cloud 
size R c ^100 pc. We assume the boundary layer containing one 
layer of similar clouds, i.e. its height is characterized by the size 
of the clouds R c . Consider an X-ray photon emitted from the 
jet with an angle 8 with respect to the jet axis in the observer 

o 

frame; for FSRQ, 8 is small, < 10 . The average number 
of the clouds that the photon is expected to encounter before 
it escapes out of this 'cocoon' is Mnc = R c /(sinOd c ), where 
d c is the mean distance between two neighboring clouds. If 
the clouds are abundant enough, d c /R c < (sin0) _1 , Mnc reaches 
unity (Mnc ^ 1) and results in X-ray absorption with almost a 
full coverage as viewed from the observer. The column den- 
sity is then N u ~ nR c N enc = 3N enc 10 21 " 22 cm" 2 . Thus, for only 
one encounter on average (Mnc=l) me resulting Nu is suffi- 
cient to account for the X-ray column densities found in most 
of the excess absorption quasars. To reach Nh ~ 10 23 cm -2 in 
RXJ1028. 6-0844, an increase in Mnc U P t0 several encounters 
is required; this may be achieved by a rather small viewing an- 
gle to the jet 8, and/or an increased cloud number density (re- 
duced d c ). For such a model, the presence of X-ray absorption 
in radio quasars might be governed by the formation of such a 
boundary layer, for which cosmic evolution is needed to explain 
the difference at high and low redshifts. 



5.2. The quasar 

The 2-50keV luminosity of 2.6 10 47 ergs _1 , after correc- 
tion for the total absorption, makes RX J 1028. 6-0844 one of 
the most luminous quasars. Furthermore, if Az < 0.1 Az©, we 
have Nu > 10 24 cm" 2 and the Thomson depth m, — 1.2cr T hMt 
reaches unity for photons with energies above a few keV (so 
that electrons bound to atoms can be treated as free); electron 
scattering becomes non-negligible and actually dominating the 
opacity for photons with E > 5 keV. An additional correction 
for such electron scattering opacity would raise the 2-50 keV 
luminosity at least a few times higher, perhaps approaching 
10 48 ergs" 1 . This will make the quasar one of the most lumi- 
nous steady objects if the radiation is isotropic. 

Fig. 5 shows the spectral energy distribution (SED) for 
RX J 1028. 6-0844, with the X-ray luminosity corrected for the 
total absorption. The UV luminosities (Luv) are somewhat un- 
certain due to the unknown dust extinction possibly imposed 
by the X-ray absorber. For a comparison, also plotted are Luv 
(open circles) for one extreme case in which the absorption is 
solely attributed to a local intervening (galactic) absorber with 
the Galactic dust-to-gas ratio; Luv then reaches 10 47 ergs _1 . 
The situation is less clear for the case of intrinsic absorption, 
as discussed above; in any case Luv must be no less than the 
observed values for which only the correction for the Galactic 
reddening was applied (filled dots). 

It has been pointed out by Fabian et al. (1999) that this object 
(and the other two z > 4 quasars) is distributed somewhat apart 
from the bulk of nearby blazars on the a rx -a ox and a ro -a ox 
planes, where a ro , a rx , and a ox are the so-called broad band 
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effective spectral indices 9 spanning between the radio, optical, 
and X-ray wavebands, respectively. For RX J 1028. 6-0844, this 
might be, at least partly, a consequence of the optical/UV ex- 
tinction imposed by the obscuring matter. The intrinsic a IO and 
a ox would flatten and steepen, respectively, and this would shift 
the object toward or into the bulk distribution of FSRQ and 
radio-selected BL Lac objects in the a ro - a rx - a ox parame- 
ter space (see Fig. 4 in Fabian et al. 1999), depending on the 
amount of dust extinction. For example, assuming local ab- 
sorption and the dust-to-gas ratio similar to the Galactic one, 
we found 10 intrinsic a ro =0.56, a rx =0.75, and a ox =l.ll. 

Both, the broad band SED and the extreme X-ray luminos- 
ity suggest that RX J 1028. 6-0844 might be a blazar type ob- 
ject; its high radio power is also of typical FSRQ. We may thus 
come to the suggestion that all these three z > 4 quasars, re- 
gardless of the strong X-ray absorption of RXJ1028. 6-0844, 
are of the same class, and the X-rays are thought to be emit- 
ted via inverse Compton radiation (Jones et al. 1974). Similar 
to GB 1428+4217 as discussed in Fabian et al. (1999), the opti- 
cal/UV luminosity of RX J 1028.6-0844 might not be dominated 
by the beamed non-thermal emission, but originate from accre- 
tion process; if this is true, the high luminosity would imply that 
black holes as massive as 10 8 — 10 9 M Q have formed within less 
than 1 billion years in the early stages of the universe. 

Owing largely to its low-energy cutoff, the X-ray spectrum 
of the quasar is consistent with that of the cosmic X-ray back- 
ground (CXB) in the 0.6-10keV ASCA band, T=1.4 power 
law (Gendreau et al. 1995) or kT =40keV bremsstrahlung (3- 
50keV, Marshall et al. 1980) with Galactic absorption (fits with 
the parameters fixed at these values resulted in reduced x 2= l l 
and 1.0 for 212 d.o.f., respectively). However, this type of ob- 
jects appears simply too rare to make significant contribution to 
the CXB. 

Finally, it is noted that the z=4.3 quasar 1508+5714 revealed 
a high rest-frame temperature 93(^4) keV (90% errors, Moran 
& Helfand 1997) when fitted with a redshifted bremsstrahlung 
with Galactic absorption model. This value is similar to that 
found for RX J 1028 .6-0844 of 132(^)keV, which apparently 
resulted from the low-energy cutoff in its spectrum. Such a 
comparison suggests possible excess absorption in 1508+5714 
too, of which the ASCA observation might not be deep enough 
to allow a detection. The 90% upper limit on the excess ab- 
sorption Nu at the quasar redshift was given as 1.3 10 22 cm" 2 
(Moran and Helfand 1997). 

6. SUMMARY 



We present the broad band 3-50 keV (source rest frame) X- 
ray observation of the quasar RX J 1028. 6-0844 at z = 4.28. The 
spectrum, with its unprecedented high quality among observa- 
tions of objects at similar redshifts, reveals an evident, substan- 
tial low-energy cutoff in excess of that due to Galactic absorp- 
tion. The spectrum is best modeled by a power law with excess 
photoelectric absorption model, yet a broken power law with 
Galactic absorption cannot be ruled out. We suggest the ex- 
cess absorption to be the most likely explanation, considering 
the spectral data at the low-energy end. The column density of 
the absorber, depending on its redshift and metallicity, ranges 
from 2.5 10 21 cm" 2 for local absorption up to 2. 1 10 23 cm -2 (so- 
lar metallicity) or 1 .6 10 24 cm" 2 (10% solar metallicity) for ab- 
sorption at the quasar redshift. It has to be noted there may 
be a systematic uncertainty in the determination of the column 
density, which tends to overestimate the absorption; however, 
we suggest this effect to be small and negligible in comparison 
to the statistical uncertainties. We attribute a possible part of 
the absorption opacity to a putative foreground spiral galaxy; 
the remaining absorption, most likely the majority part, seems 
not to be easily explained by cosmological intervening systems. 
Future observations with XMM and Chandra would be able to 
measure the redshift and other parameters of the absorber by 
detecting possible iron absorption edges above 7 keV. 

The unabsorbed quasar continuum is well described by a 
simple power law with T = 1.67^^ extending up to 50keV 
in the source rest frame. The absorption corrected luminos- 
ity reaches as high as 2.6 10 47 ergs _1 , regardless of possible 
electron scattering effect. No statistically significant X-ray flux 
variations were found during the relatively short observational 
interval, neither by a comparison with the early ROS AT Survey 
observation. Considering its extreme apparent luminosities in 
the X-ray and radio bands, as well as the broad band SED, we 
suggest RXJ1028. 6-0844 to be a blazar type object, similar to 
the other two z > 4 quasars observed with broad band X-ray 
spectroscopy. 
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Table 1 

asca observation and data screening 



Item 


SIS0 


SIS1 


GIS2 


GIS3 


Elevation angle 8 
C.O.R." 
Exposure (sec) 
Counts' 

Background counts' 1 
Net count rates 


15°/30° 
8 

59061 
2179 

675 (1.14) 
2.54±0.10 


15°/25° 
8 

60533 
1831 

606(1.08) 
2.02±0.09 


15° 
8 

55085 
1620 
657 (2.0) 
1.74±0.08 


15° 
8 

55090 
1807 
786 (2.0) 
1.85±0.09 



"For SIS the second number refers to the bright Earth elevation angle. 

b Magnetic cutoff rigidity in GeV/c 

'Overall counts from the source extracting region 

''Background counts normalized to the source extracting region, which are the ex- 
tracted background counts divided by the ratio of the background-to-source region 
areas (in bracket). 

'Background subtracted count rate in units of 10~ 2 ctss~' in the 0.6-10keV band 
for SIS and the 0.7-10 keV band for GIS. 
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Table 2 

Results of X-ray spectral fits 

Power law + local neutral absorption 
Detector N H a T norm. b x 2 /dof 



GIS 

SIS0+SIS1 

SIS+GIS 

GIS 

SIS0+SIS1 
SIS+GIS 



0.0459 (fix) 

0.0459 (fix) 

0.0459 (fix) 

30 +012 
u - JU -o.io 

26 +006 
30 +0 05 



1.51 ±0.06 
1.24 ±0.04 
1.34 ±0.03 
1 76 +012 

i - / D -0.07 

1.56 ±0.09 
1 67 +007 

1 - o/ -0.04 



1.98±0.12 
1.38 ±0.06 
1.69 ±0.07 

9 C 1 +0.51 

2 06 +0 ' 24 
2 61 +0 ' 24 



110/100 
102/112 
221/211 
104/99 
84/1 1 1 
187/210 



Power law + excess absorption at z = 

A Z /A ZQ N£ c r 


4.28 andJVg al (SIS+GIS) 
norm. x 2 /dof 


1.0 (fix) 
0.1 (fix) 


21 1 +5 - 3 1 72+ 010 

Z1 -M.2 i -' z -0.07 
1^+31 i 79+0.04 


2.86!g^ 187/210 
2.81 ±0.50 186/210 


Bremsstrahlung + local neutral absorption (SIS+GIS) 
redshift N H T (keV) norm. c x 2 /dof 


4.276 (fix) 
4.276 (fix) 


0.0459 (fix) 13 

n?1 +0.04 ca 7+12.2' 
u - zi -0.03 J ^-'-9.2 


3.53 ±0.16 210/211 
3.75^3;^ 187/210 


Broken power law + local neutral absorption (SIS+GIS) 

£break d N H F lov/ Thigh X 2 /dof 


1 7 +0 ' 2 

2 +mf 


0.0459 (fix) 0.70^ 
0.28 ±0.07 1.58^| 


1.58!g?f 190/209 
1.67!°;^ 187/208 



a Column density of hydrogen in units of 10 cm 

b Nomalization in units of 10~ 4 photons s" 1 cm" 2 keV" 1 at 1 keV. For 
joint SIS+GIS fits, the values given are of GIS, and those of SIS are 
less by within 10% in general. 

c Nomalization as given in XSPEC in units of 
10" 14 /(4ttD 2 ) J n^n^dV, where D is the luminosity distance to 
the source (cm) and n e and «h the electron and hydrogen densities 
(cm -3 ), respectively. 

d The break energy of broken power law in keV. Un-specified un- 
certainties mean the value is unconstrained. 
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0.5 1 2 5 

Energy (keV) 



FIG. 1 . — Upper panel: Joint spectral tit of a power law with fixed Galactic absorption to the SISO (circles), SIS1 (squares), and the combined GIS2+GIS3 (stars) 
data in the restricted 2.5-9.5 keV energy band (with the normalizations untied). The model values (convolved with the detector responses) are represented by lines 
from upper to lower (at low energies) for SISO, SIS1, and GIS2+3, respectively. Lower panel: Residuals showing deviations of the data from the extrapolation of 
the fitted model down to the low energy range for all the detectors, respectively. 




0.5 1 2 5 

Energy (keV) 

FIG. 2. — Residuals for the best-fit model — a power law with free absorption column density to the joint SIS and GIS data. The SIS data below 0.7 keV are not 
used in the fitting but plotted for a comparison. The plot symbols are the same as in Fig.l. 



X-ray absorption in the high-z quasar RXJ1028. 6-0844 
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0.1 1 10 

N„ (10 Z2 cm -2 ) 

FIG. 3. — Confidence contours of fitted absorption column density and photon index for the total absorption at redshift z=0 and for the excess absorption (with 
additional fixed Galactic absorption) assumed at the quasar redshift z=4.28, respectively. The contours are at the 68%, 90%, and 99% confidence, respectively, for 
two interesting parameters. The elemental abundances are assumed to be the solar values. The best-fit values are indicated by crosses. The Galactic column density 
and a conservative 30% uncertainty region are indicated by dashed and dotted lines, respectively. 




°0 1 2 3 4 



Redshift 

FIG. 4. — Contours of allowed excess absorption A^ xc and the redshift of the absorber for two interesting parameters at the 68%, 90%, and 99% confidence levels, 
respectively. The elemental abundances are assumed to be the solar values. 
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FIG. 5. — Spectral energy distribution of the quasar in the source rest frame, in which the luminosities are at the emitting frequencies. The X-ray luminosities 
are corrected for the total absorption. The RASS measurement (triangle) is converted from the count rate assuming the spectral model as the extrapolation of the 
best-fit absorbed power law derived in this work, and at a frequency corresponding to 1 keV in observer frame. The UV luminosities are obtained from the optical B- 
and R-band (non-simultaneously) photometry data (Zickgraf et al. 1997), with the correction for dust extinction applied for: (a) the Galactic extinction only (filled 
dots); (b) the Galactic plus an addition column of a local X-ray absorber (z < 0.1) assuming the Galactic dust-to-gas ratio (open circles). In the case of the excess 
absorption occurring at a high redshift or intrinsic to the quasar, the intrinsic UV luminosities are uncertain due to the unknown extinction; however, they must be 
no less than the observed values for which only the correction for the Galactic reddening was applied (filled dots). 



